InN dielectric function from the midinfrared to the visible range 
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The dispersion of the dielectric function for wurtzite InN is analytically evaluated in the region 
near the fundamental energy gap. The real part of the dielectric function has a logarithmic singular- 
ity at the absorption edge. This results in the large contribution into the optical dielectric constant. 
For samples with degenerate carriers, the real part of the dielectric function is divergent at the 
absorption edge. The divergence is smeared with temperatures or relaxation rate. The imaginary 
part of the dielectric function has a plateau far away from the absorption onset. 
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Recently, InN has attracted considerable attention due 
to its potential application as other Ill-nitrides, but es- 
pecially owing to the small energy band gap 2e g of about 
0.7 eV observed [H, [2, & El, HI m contrast to the value of 
1.9 eV established for the last 20 years. The small band 
gap value corresponds with a small effective electron mass 
to* » 0,07m 0,0,11. 

For future progress in the research field, reliable mate- 
rial parameters are derived from the most widespread ah 
initio electronic-structure calculations. However, these 
methods do not present analytical results and lead some- 
times to contradictions Q , whether the 4d bands are in- 
cluded in the core or are not. Therefore, the k • p Hamil- 
tonian is used to clarify the physical content. In the 
corresponding Kane model for the wurtzite case, the 
conduction-band and the valence-band are constructed 
from the \s) and \y), and \z) states at the T— point. 

In Fig. QJ,, the scheme of the valence-band splitting in 
InN is shown under the crystal field Acr and the spin- 
orbit interaction Ago- According to experimental data 
fli"| and calculations [|[ , this splitting has a value on the 
order of 0.02-^0.06 eV, i. e. it is small in comparison with 
the band gap, and can be ignored in calculations of the 
integral properties as the optical absorption. Therefore, 
the Kane model can be substantially simplified while us- 
ing in calculations of the dielectric function (DF). 

In this paper with the help of the simplified Kane 
model, we evaluate analytically the DF for wurtzite InN 
in the range 10 ~ 0.3 -j- 4 eV, where the absorption is 
dominated by optical transitions from three highest va- 
lence bands into the lowest conduction band. As well 
known, the imaginary part of the DF has the square- 
root behavior near the absorption edge in the case, if 
the conduction-band is empty (i), and the step-like be- 
havior in the case, if the carriers in the conduction band 
are degenerate (ii). In the paper [l2j|, the ab initio cal- 
culations of the imaginary part were presented for the 
InN polymorphs (with wurtzite, zinc-blend, and rocksalt 
structures) in the case (i), whereas the real part was re- 
stored using the Kramers-Kronig relations. In the paper 
[l2| . the imaginary part was also estimated within the 
Kane model, but these results are misleading, because the 
optical-transition matrix elements are incorrectly cvalu- 



FIG. 1: (a) scheme of the valence-band splitting under the 
crystal field and the SO interactions; (b) the electron band, 
the heavy-hole and light-hole bands near the T— point. 



ated. 

We find, that the real part of the DF has at the ab- 
sorption edge the kink-like singularity in the case (i) and 
the logarithmic divergence in the case (ii). Such singu- 
larities have been previously obtained for graphene [13j 
and for IV- VI semiconductors [HI]. The singularities are 
smeared with temperatures or carrier relaxation. The 
excitonic effects was not observed in InN, since they are 
suppressed by the carrier relaxation (see 0, [l2|, [l5[ for 
their estimation). 

The effective Hamiltonian of the simplified anisotropic 
Kane model is given as the matrix 4x4: 
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where P\ ^ P 2 are the interband momentum-matrix el- 
ements with the velocity dimension (we put h = 1 in 
the intermediate formulas). Quadratic terms in the mo- 
mentum k can be written [16j in the main diagonal of 
the matrix JT]), as well as in the terms, connecting the 



states at the valence band top, — e g . We omit them, 
because their contribution in the DF is on the order of 
1/ln (e a t/E g ) <C 1, where e a t is the energy of the atomic 
scale. 

This Hamiltonian ([I]) gives rise to the eigenvalues 



ei,4 = ±[el + Pfkl + P 2 k 2 ± } 1/2 , k\ = k'l 



(2) 



corresponding with the conduction band and the light- 
hole band, and the twofold eigenvalue 



-2,3 



(3) 



for the heavy- hole band (see Fig. [T]d). 

The effective masses at the conduction-band bottom 



arc 



m,| = e g /P£ , m± = e g /P 2 (4) 

for the longitudinal and transverse directions respectively 
to the main z— axis. Comparison with the values my = 
0.065m , m± = 0.068m , and 2e g = 0.69 eV obtained 



in the experimental data analysis gives Pi 
cm/sec, P2 = 0.945 x 10 s cm/sec. 
The velocity operator has the form 



0.966 x 10 s 



(0 iP 2 e x iP 2 e y iPie z 
-iP 2 e x 
-iP 2 e v 
-iPie z 

with orthonormal vectors chosen along the coordinate 
axes. 

Let us define the velocity matrix in the representation 
diagonalizing the Hamiltonian JT]): 

v = U^vU. 

Using the eigenfunctions of the Hamiltonian, one finds 
the matrix U : 



U 



I — Ol d 2 

ipx/ni ~Py/P± PxPz/PxP iPx/n 2 

ipy/nt p x /p± P y Pz/p±P ip y /n 2 

\ iPz/ni ~P±/p ip z /n 2 



where we use the notations p z = P\k z ,p x . y = 
P 2 k x , y ,p = \/y z + p± , ni,2 = y/2ei(ei ± e g ) , 01,2 = 
\J (ei ± £g)/2e\, and £\ is given in Eq. (fj)). 

We note that the diagonal velocity-matrix elements in 
this representation coincide with the derivative of the 
eigenvalues 

de l 

(«M.) = ^ 
and we write the off-diagonal elements 



(l|v|3) = 
(l|v|4) = 

(3|v|4) = 



-1P1P2 



£l\/£l-£g 

-iPiP 2 



k 



y / 2e 1 (e 1 + e g 



k 



~{kx&X ~t~ kyGy^) kj_G 2 



These matrix elements enter the general quantum- 
mechanic formula for the dynamic conductivity o a 0{oj) 
derived in the paper Then, we obtain the DF with 
the help of the relation 



3/3(0;) = 1 + 4iri<7 al3 (uj)/uj . 



(6) 



Due to the symmetry, the off-diagonal tensor components 
of the DF vanish and there are only two independent 
components e zz (uj) and e xx (uj) — e yy (tu). 

The DF is separated into the intraband and intcrband 
parts. The intraband term contains only the diagonal 
velocity-matrix elements and has the Drude-Boltzmann 
form. For instance, we obtain for the degenerate elec- 
trons: 



(7) 



where the chemical potential ^, the relaxation rate v and 
the photon frequency to are written in the common units. 

Neglecting the carrier relaxation, we can write the in- 
tcrband term for the extraordinary component of the DF 
in the form 



2e 2 

"H = l + -2" 
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(8) 



[/(-ei) - f(-e s 
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(£l-£ 5 )[(£l-£ ff ) 2 -(w + ^) 2 ] 



where /(e) = l/(cxp [(e — p)/T] — 1) is the Fermi func- 
tion. For the pristine semiconductor at low temperature, 
the conduction band is empty, but the chemical potential 
fx can be higher than the conduction band bottom e g in 
the case of doping (see Fig. [TJa). 

The different terms in the braces present the verious 
optical transitions: first, between the light-hole band and 
the conduction band, second, between the heavy-hole 
band and the conduction band, and third, between the 
light- and heavy-hole bands. The infinitesimal 8 in the 
denominators of Eq. ([5]) defines the bypass around the 
poles. These bypasses give the imaginary part of the DF, 
whereas the principal values of the integrals yield the real 
part. 

Transforming the integration variables 
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where 



sin 0,(9) 



Pi d(e u 8,<p) 
we integrate over the angles and 
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/(- £ ) - /(-e fl ) 



(e + e 3 ) 2 - (w + z<5) 2 (e-e 3 ) 2 -(w + ^) : 



The integral presenting the real part of the DF diverges 
logarithmically at the upper limit. Since the leading con- 
tribution arises from the values e ~ (n,ui), the integral 
can be cut off at the atomic value of energy e at , where 
our k • p expansion becomes incorrect. The imaginary 
part is easily evaluated for zero temperatures. For in- 
stance, we find for the case \i > s g , when electrons fill 
the conduction band, 



Ime» ter H 



4e 2 F! e g 
3hP2 I 2uj 2 
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where the step 9— function conveys the condition for the 
interband electron absorption. 

Let us emphasize, that the band edge for the optical 
transitions into the conduction band from the light-hole 
band at ui = 2\i is higher than the edge for the tran- 
sition from the heavy-hole band at u> = e g + (i. With 
increasing the free electron concentration, both edges 
demonstrate the blue Burstein-Moss shift. At zero tem- 
peratures, the chemical potential /i, measured from the 
midgap is determined by the free-electron concentration: 
n = (M 2 -4) 3/2 /3^ 3 PiP 2 2 . i 

If the electrons are absent in the conduction band, 
— e g < [i < e g , the imaginary part of the DF is given 
in Eq. (|10p with substitution /i — > E g . Far away from the 
absorption edges, where u> 3> E g + fi, the imaginary part 
demonstrates the plateau-like character with 



max Im e l 



4e 2 Pi 

3hpJ 



(11) 



The plateau noticed also in the paper [l2| and for the 
A4B6 semiconductors in [l4| is a consequence of the lin- 
earity of the electron dispersion at the energy larger in 
comparison with the energy gap. 

The real part of the DF contains the following contri- 
butions. The transitions between the heavy-hole bands 
and the conduction band give 
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F{ x ) = — — In - V 
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if x = Eg/uj < 1/2 and 



(13) 



+ 1. 



F(x) = —\/2x — 1 arctan ■ 



}y/2x 



Ll + £g 



M 2 ~ £ 2 Q 



if x > 1/2. 

The transitions between the light-hole band and the 
conduction band contribute 
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for x < 1/2 and 
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for x > 1/2. 

We find that the real part of the DF as a function of 
ui takes at x — 1/2 the maximal value for fj, = e g : 
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If the carriers appear in the conduction band, fi > s g , 
the real part of the DF has a logarithmic singularity 
smeared with temperature or carrier relaxation. For 
small relaxation rate v in comparison with the photon 
frequency u>, we have in Eqs. (fT0|) . (fT3|) . (fT4)l to substi- 
tute 



- H — arctan[(cj — cj a t)/2f] (17) 

2 7T 



(<J - U> at ) 

(uj-uj at ) 2 (uj-uj at ) 2 + (2v) 2 



where ui a t is the absorption edge equal to s g + fi or 2fi 
for the corresponding transitions. If temperature plays a 
more important role, we should put T instead of v in Eq. 

(HZD- 

So far the extraordinary component t zz was presented. 
The ordinary component e xx differs only in the factor 
P2/P1, which equals 0.98 for the experimental values of 
the effective masses (fj|. 
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FIG. 2: The real and imaginary parts of the DF versus the 
photon frequency (in units of the gap 2e g — 0.69 eV) for 
various free-electron concentrations, corresponding values of 
the chemical potential /x are 1.01, 1.57, 2.79 (in units of 2e 9 = 
0.69 eV); relaxation rate v — 0,01 /j,. 



Now the value of the cutoff parameter e a t is only 
needed to calculate the DF. To estimate this value, we 
can use the energy arising in the Kane model while the 
quadratic terms are taken into account. According to 
estimations [1, [l2[, this energy ranges from 8 to 15 eV. 
We take the intermediate value e at — 10 eV plotting Fig. 
[2l where our theoretical results are shown. The maxi- 
mum value of the real part, Eq. (|16[) . is found to equal 
6.91 and the imaginary part of the DF takes the value 



3.16 on the plateau, Eq. (fTT|) , The corresponding values, 
2.5 -7- 3.5, obtained from experiments [13, El, [3 and cal- 
culated from the first principles [12] are on the plateau in 
the frequency range 1.5 -j- 4.0 eV. That agrees very well 
with Fig. [2] (right panel) . The experiment [H[ finds the 
value about ~ 9 for the maximum of the real part. The 
estimation [l8| of the dielectric constant gives = 6.7. 
The ab initio calculations [l2l , |20| find correspondingly 
in these two papers e^xx) = 7.03 and 7.16, as well as 
e oa (zz) = 7.41 and 7.27. The agreement with our Fig. 
[2] (left panel) is excellent again. In our calculations, the 
maximum of the real part for the large carrier concen- 
tration increases logarithmically with decreasing of the 
relaxation rate. Plotting Fig. [2j we take v = 0.01/j, for 
various carrier concentrations. 

Concerning the dielectric constant t m , we keep in mind 
that the phonons contribute into its value. This contri- 
bution can be estimated as lo^q/uj 2 , where u>to is the 
transverse phonon frequency which is much less than the 
photon frequency considered here. Therefore, the phonon 
contribution into eoo should be considered as negligible. 

In conclusions, we find analytically that the real part 
of the DF contains a singular contribution from the inter- 
band optical transitions. It presents the large logarithmic 
term to the optical dielectric constant. While increasing 
the frequency, we obtain the dispersion of the dielectric 
function. Near the edge of the intcrband absorption, a 
peak appears in the real part of the DF for degenerate 
electrons filled the conduction band if the relaxation rate 
is large enough. 
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